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1.  INTRODUCTION 


The  instability  and  eventual  breakup  of  shaped-charge  jets  have  been  previously  observed 
in  many  experimental  studies.  These  jets  are  commonly  formed  from  the  collapse  of  thin 
metallic  conical  liners,  when  exposed  to  the  detonation  of  a  high  explosive.  After  formation, 
shaped-charge  jets  stretch  with  a  reasonably  uniform  velocity  gradient,  when  measured  along 
their  length.  Apparently,  this  stretching  motion  is  unstable,  because  after  a  specified  time 
period  (the  breakup  time)  the  jet  will  neck  down  in  places  along  its  surface  and  break  into 
many  fragments. 

Interest  in  the  breakup  phenomenon  observed  in  shaped-charge  jets  has  prompted  a  num 
ber  of  theoretical  investigations  into  the  stability  of  elongating  jets.  Recent  two-dimensional 
calculations  by  Curtis  (1987),  Pack  (1988),  and  Romero  (1989),  in  which  the  jet  was  as¬ 
sumed  to  be  perfectly  plastic  and  obey  the  Von-Mises  yield  criterion,  have  been  successful  in 
reproducing  some  of  the  behavior  observed  in  experiments.  These  studies  also  accounted  for 
inertial  effects,  since  they  are  significant  for  shaped-charged  jets,  where  strain  rates  as  high 
as  10s  sec-1  are  commonly  observed.  The  results  indicated  that  uniformly  elongating  plastic 
jets  are  almost  always  unstable  when  a  perturbation  is  introduced.  Romero  concluded  that 
the  most  unstable  disturbances  for  jets  with  large  strain  rates  have  initial  wavelengths  that 
ar°  on  the  order  of  the  jet  diameter.  This  prediction  is  consistent  with  the  experimental  ev¬ 
idence  for  shaped-charge  jets,  where  roughly  equal-sized  fragments  have  been  observed  after 
the  jet  breakup  has  initiated.  A  later  study  by  Littlefield  and  Powell  (1990)  which  also  in¬ 
cluded  electromagnetic  effects,  discussed  in  some  detail  the  physical  mechanisms  responsible 
for  this  breakup. 

Recent  attention  has  focused  on  the  development  of  methods  for  enhancement  of  the 
instabilities  in  shaped-charge  jets,  in  order  to  reduce  the  breakup  time.  One  proposed 
method  involves  the  introduction  of  an  axial  electric  current  in  the  jet.  By  passing  an  axial 
current  through  the  jet,  it  is  hoped  that  instabilities  similar  to  those  found  in  plasma  columns 
(Kruskal  and  Schwartzschild  1954;  Tayler  1956)  might  be  activated.  Recent  experiments 
have  demonstrated  that  electrical  currents  may  have  an  effect  on  the  breakup  time  of  jets 
Moreover,  the  study  by  Littlefield  and  Powell,  in  which  a  linear  perturbation  analysis  was 
employed  to  examine  the  stability  of  infinitely  conducting  jets,  indicated  that  electrical 
currents  may  shorten  the  breakup  time  of  the  jet.  Further  study  is  needed,  however,  before 
any  definite  conclusions  can  be  drawn.  The  applicability  of  the  analysis  by  Littlefield  and 
Powell,  for  example,  is  limited  due  to  the  assumption  of  infinite  conductivity.  A  typical 
copper  shaped-charge  jet  (Powell  and  Littlefield  1990)  has  a  conductivity  of  about  107  (ohm- 
m)_l.  The  magnetic  Reynolds  number  associated  with  the  jet,  therefore,  is  about,  ten,  which 
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is  not  large  enough  to  justify  the  infinite  conductivity  assumption. 

The  present  analysis  extends  the  previous  work  by  Littlefield  and  Powell  to  include  jets  of 
finite  conductivity.  As  in  previous  analyses,  the  jet  is  assumed  to  be  infinitely  long,  uniformly 
elongating  in  the  axial  direction,  and  perfectly  plastic.  Heat  conduction,  dissipation,  and 
Joule  heating  effects  are  also  neglected,  which  permits  the  omission  of  the  conservation  of 
thermal  energy  and  simplifies  the  analysis.  In  practice,  this  assumption  restricts  the  electrical 
current  to  levels  that  are  not  laige  enough  to  cause  large  scale  changes  in  the  thermal  energy. 
Two  different  base  solutions  for  the  velocity,  pressure  and  magnetic  field  distributions  are 
determined  for  the  jet,  as  well  as  the  surrounding  vacuum,  under  axisymmetric  conditions. 
In  the  first  solution,  an  axial  electric  current,  assumed  to  be  initially  on  the  surface,  is 
permitted  to  diffuse  into  the  jet  interior  over  time.  The  second  solution,  on  the  other  hand, 
assumes  an  initial  body  current,  so  that  the  axial  current  density  is  constant  throughout 
the  cross  section.  While  the  first  solution  is  characteristic  of  the  actual  current  distribution 
in  jets,  the  second  solution  provides  useful  insights  about  the  stability  and  facilitates  the 
discussion.  Small  axisymmetric  disturbances  of  these  base  solutions  are  then  considered,  and 
the  equations  governing  their  time  evolution  are  derived  using  linear  perturbation  theoiy. 
Solutions  to  these  equations  are  used  to  evaluate  the  stability  of  the  jet. 

The  arrangement  of  this  paper  is  as  follows.  In  Sec.  2,  the  physical  model  is  described 
and  the  governing  equations  in  the  jet,  as  well  as  the  vacuum,  are  presented.  Base  solutions 
to  these  equations  are  given  in  Sec.  3  for  the  different  initial  current  distributions.  The 
first-order  equations  corresponding  to  these  base  solutions  are  derived  in  Sec.  4  using  ljneai 
perturbation  theory.  In  Sec.  5,  these  equations  are  solved  and  the  stability  characteristics  of 
the  jet  are  discussed.  Finally,  Sec.  6,  contains  some  general  conclusions  and  observations. 

2.  MATHEMATICAL  MODEL 

Consider  an  axisymmetric  shaped-charge  jet  as  depicted  in  Fig.  1.  The  jet  is  assumed  to 
be  uniformly  elongating  and  perfectly  plastic.  A  constant  axial  electric  current  I  is  pass(  d 
through  the  jet  initially  and  permitted  to  diffuse  with  time.  The  radius  at  any  point  alon' 
the  jet  boundary  is  denoted  by  r(,(z,<),  and  the  outward  unit  normal  vector  at  this  point 
by  n.  The  governing  equations  of  motion  in  the  jet,  which  may  be  found  elsewhere  (Poweli 
and  Littlefield  1990;  Littlefield  and  Powell  1990)  are  repeated  here  for  convenience.  Let  the 
velocity  be  given  by  V,  the  pressure  by  p ,  the  deviatoric  stress  tensor  by  5,  the  current 
density  by  J,  the  electric  field  by  E,  the  magnetic-induction  field  by  B,  the  density  by  p ,  the 
conductivity  by  a,  and  the  magnetic  permeability  by  po-  Then  the  governing  equations  are 
given  by 
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VV  =  0, 


(1) 


fdv 

p[-dr 

+  V-VvJ  =  -Vp+V-3  +  JxB, 

(2) 

E  +  V  xB  =  J/<r, 

(3) 

VB  =  0, 

(4) 

V  x  B  =  hqJ, 

(5) 

and 


V 


(6) 


Figure  1.  Model  for  jet  stability  calculations. 

Equations  (1)  and  (2)  are  the  incompressible  continuity  and  momentum  equations,  re¬ 
spectively;  Eq.  (3)  is  Ohm’s  Law,  and  Eqs.  (4)  -  (6)  are  Maxwell’s  equations. 
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The  present  analysis  assumes  the  idealized  motion  of  the  jet,  as  well  as  disturbances  to 
this  motion,  are  axisymmetric.  Furthermore,  Jg  and  Eg  are  assumed  zero.*  Consequently, 
the  component  forms  of  the  governing  equations  may  be  simplified.  The  continuity  equation, 
tor  example,  reduces  to 


Id  dv2 

A(r’’)  +  ¥  =  0' 


(7) 


where  vT  and  vz  are  the  radial  and  axial  components  of  the  velocity,  respectively.  Ampere’s 
law,  given  in  Eq.  (5),  yields  the  two  scalar  relations 


and 


;|(r  B)  =  /W„  (!>) 

where  JT  and  Jz  are  the  radial  and  axial  components  of  the  current  density,  respectively, 
and  B  is  the  azimuthal  component  of  the  magnetic  field.  With  use  of  Eqs.  (8)  and  (9),  the 
component  forms  of  the  momentum  equation  become 


dvr  dvr 

+  IV—  +  Vz  — 
Or  dz 


dp  dsTT  sTT  sgg  dsrz  1  B  d 
dr  dr  r  dz  r  po  dr 


(10) 


and 


P 


( 


dvz  dvz 

W  +  v^  +  v 


dp  dsrz  srz  dsTZ  B  dB 

di +  IT  +  T  + 17  ”  7o  a7’ 


(ii. 


where  sap  denotes  the  corresponding  component  of  Ohm’s  Law  yields  the  two  component 
equations 


a(ET  —  vzB)  =  Jr 


(12) 


and 

"The  governing  equations  allow  for  nonzero  Jg  and  Eg,  but  it  may  be  shown  these  quantities  tend  to  zero  as  time  progresses 
anyway.  Consequently,  the  contribution  of  nonzero  Jg  and  Eg  to  the  stability  characteristics  of  the  jet  is  believed  to  be  small. 
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o(Ez  +  vrB)  =  Jz, 


(13) 


where  Er  and  Ez  are  the  radial  and  axial  components  of  the  electric  field,  respectively.  Eq.  (4) 
yields  information  uncoupled  from  the  present  analysis,  but  (6)  yields  the  single  equation 

3EL_dE1  =  _dB_ 

dz  dr  dz 

The  scalar  forms  of  Eqs.  (1)  -  (6),  given  in  Eqs.  (7)  -  (14),  still  have  more  unknowns 
than  equations.  A  constitutive  relationship  must  also  be  specified,  relating  the  deviatoric 
stress  to  material  properties  and  other  variables  in  the  jet.  In  the  present  analysis,  the  jet 
is  assumed  perfectly  plastic.  This  requires  the  deviatoric  stress  to  be  proportional  to  tin 
strain  rate,  which  under  axisymmetric  conditions  yields 

srr  =  2A^,  (15) 


(16) 


Szt 


(H) 


and 


=  A 


fdv^  dvz\ 
\dz  +  dr  )' 


(18) 


where  the  coefficient  of  proportionality  A  is  formulated  to  satisfy  the  Von  Mises  yield  crite¬ 
rion,  given  by 


(Srr  -  seg )2  +  (see  -  SZZ )2  +  (szz  -  srT)2  +  64  =  2F2,  (19) 

and  Y  is  the  constant  yield  strength.  Solving  for  A  from  Eq.  (19)  with  use  of  Eqs.  (15)  -  (18) 
gives 
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Application  of  these  constitutive  relations  to  the  momentum  equations  given  in  Eqs.  (10) 
and  (11)  yields 


and 


(dv2  dv2  dv2 ' 

(ar  +  'var  +  ,’-a7 


dp  «  / d2vz  1  dvz  d2vz\ 
dz  \  dr 2  r  dr  dz 2  J 


fdvT  du2\  B  dB 
dz  dz  ^  dr  \  dz  +  dr  )  p0  dz  ' 


(22) 


with  A  given  by  Eq.  (20).  It  is  also  convenient  to  combine  Eqs. (8)  -  (9)  and  (12)  -  (14)  into 
a  single  equation  for  B,  which  yields 


dB  ^  dB  vtB  dB  1  fd2B  1  dB  B  d2B\ 

dt  +Vr  dr  r  +  Vz  dz  ~  an 0  (  <9r2  +  r  dr  r2  +  dz2  )  '  ^ 

Eqs.  (7)  and  (20)  -  (23)  represent  five  equations  for  the  five  scalar  unknowns  (vr,  v2 ,  p,  A, 
and  B)  which  govern  the  motion  in  the  interior  of  the  jet. 

The  relevant  conservation  equations  also  must  be  satisfied  in  the  vacuum  surrounding 
the  jet.  In  the  vacuum,  Eq.  (1)  is  meaningless  and  Eq.  (2)  is  identically  satisfied  since  />. 
p,  and  J  are  zero.  Eqs.  (3)  -  (6)  may  be  simplified  by  setting  J  =  cr  =  0  and  requiring 
axisymmetry.  The  resulting  equation  for  Bv,  the  0  component  of  B  in  the  vacuum,  is  given 
by 


d2Bv  1  dBv 

dr2  ^  r  dr 


Bv  d2Bv 
7*  +  dz2  ~  °' 


(24) 


The  remaining  variables  in  the  vacuum  are  either  zero  or  may  be  determined  by  back  sub¬ 
stitution. 

The  application  of  these  equations  implies  that  certain  conditions  must  be  satisfied  at 
the  interface  between  the  vacuum  and  the  jet.  In  particular,  mass  and  momentum  must  be 
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conserved  and  Maxwell’s  equations  must  be  satisfied  across  the  interface.  These  conditions 
may  be  derived  systematically  (Kruskal  and  Schwartzschild  1954)  by  integrating  the  rele¬ 
vant  equations  across  the  boundary,  of  presumed  thickness  S,  in  a  direction  normal  to  the 
boundary,  then  taking  the  limit  as  6— >0.  Let  dl  be  an  element  of  arc  length  normal  to  the 
surface.  Then,  the  pertinent  equations  are  multiplied  through  by  dl  =  nTdr  +  n2dz,  where 
nr  and  nz  are  the  r  and  z  components  of  n,  and  integrated  across  the  boundary.  Integration 
of  the  quantity  |£,  for  example,  produces  the  term  nr  <  F  >,  where  the  brackets  denote 
the  change  in  the  quantity  across  the  boundary.  Similarly,  the  term  produces  n2  <  F  >. 
Finally,  the  main  contribution  of  the  partial  time  derivative  of  a  function  F  arises  from  the 
motion  of  the  boundary,  so  that  integration  of  the  term  produces  —  v  •  n  <  F  >,  where 
v  is  the  velocity  of  the  interface. 

Application  of  the  integration  procedure  to  Eq.  (1)  conserves  mass  and  gives 


drb  drb 

ar  +  tv97  =  lv 


(2.5) 


when  evaluated  at  r  =  rb.  The  integration  of  Eq.  (2)  conserves  momentum  at  the  interface 
and  yields 


at  r  =  r( ,,  where  Eqs.  (15)  -  (18)  have  been  substituted  for  the  components  of  3.  Here  the 
radial  and  axial  components  of  the  unit  normal  vector  n,  given  by  nT  and  nz,  are  easily 
defined  in  terms  of  the  equation  of  the  surface,  r  =  rb.  This  gives 


_  V(r  -  rb) 
~  |V(r-r*)| 


1  + 


(28) 


where  er  and  e,  are  the  unit  normal  vectors  in  the  r  and  z  directions,  respectively.  Similarly, 
the  integrations  of  Eqs.  (3),  (4),  and  (6)  may  be  performed,  but  yield  either  an  identity 
or  information  which  is  not  required  for  the  analysis.  The  integration  of  Eq.  (5),  however, 
requires  B  to  be  continuous  at  the  interface,  or 
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Bv  -  B  =  0 


(29) 


at  r  =  rj.  Eqs.  (25)  -  (27)  and  (29),  together  with  the  requirement  that  vT,  v2,  p  and 
B  all  remain  finite  at  the  centerline,  define  a  complete  set  of  boundary  conditions.  If  the 
initial  conditions  are  specified,  these  boundary  conditions  are  sufficient  to  solve  the  governing 
equations  in  the  jet  and  the  vacuum. 


3.  BASE  SOLUTIONS 


The  unperturbed  motion  of  the  jet  is  characterized  by  a  uniform  velocity  gradient  in  the 
axial  direction.  This  velocity  gradient  determines  both  the  axial  and  radial  local  velocities, 
as  well  as  the  position  of  the  interface  at  any  given  time.  In  addition,  a  constant  axial  electric 
current  I  is  passed  through  the  jet  and  permitted  to  diffuse  with  time.  This  idealized  current 
distribution  induces  azimuthal  magnetic  fields  both  in  the  jet  and  in  the  surrounding  vacuum, 
which  interact  to  create  Lorentz  forces.  These  electromagnetic  forces,  as  well  as  inertial  forces 
which  arise  from  the  velocity  gradients,  affect  the  idealized  stress  distribution  in  the  jet. 

The  unperturbed  velocities  and  radius  of  the  jet  have  been  presented  many  times  (Powell 
and  Littlefield  1990;  Curtis  1987;  Pack  1988;  Romero  1989;  Littlefield  and  Powell  1990),  but 
will  be  repeated  here  for  convenience.  The  subscript  zero  will  be  used  to  denote  these 
idealized  quantities.  The  axial  and  radial  velocities  are  given  by 


and 


(30) 


(31) 


where  t)  is  the  initial  strain  rate  and  r  =  1  +  pt.  This  velocity  distribution  requires  the 
unperturbed  value  of  A  to  be  A0  =  Yr/3r].  The  radius  at  the  boundary  is  given  by 


(32) 


where  a  is  the  initial  radius  of  the  jet.  The  magnetic  field  and  stress  profiles  depend  on  these 
velocities,  and  on  the  current  distribution  in  the  jet.  The  magnetic-induction  field  in  the 
vacuum,  for  example,  is  required  by  Ampere’s  Law  to  be 
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(33) 


ov  _  HO1 
0  "  2ttt  ’ 

which  also  satisfies  Eq.  (24).  The  zero-order  pressure  distribution  in  the  jet  is  determined 
from  the  integration  of  Eq.  (21)  which,  after  substitution  of  the  velocities  given  in  Eqs.  (30) 
and  (31)  and  use  of  the  boundary  condition  in  Eq.  (26),  yields 


-(a2  -  rr2)  -  ]r  +  —  f 
3  UO  Jr 


Y  1  p*  B0(u,  t  )  3,  n(  „  , 

— +  —  /  - —  [uB0(u,t)  ]du, 

3  Uo  J T  u  Ou 


where  B0  is  the  zero-order  magnetic  field  in  the  jet. 

This  magnetic  field  must  satisfy  Eq.  (23),  which  may  be  simplified  to  yield 


dB0  rjr  dBo  t]B0  1  fd7B0  Id  B0  B0 

r\  4”  -  I  r\  n  4”  ^  0 


dt  2 r  dr  2  r  (Tfi0  \  dr2  r  dr  r7/ 

The  advection  term  in  Eq.  (35)  may  be  completely  eliminated  by  a  coordinate  transformi  - 
tion  which  deforms  in  accordance  with  the  stretching  motion  of  the  jet,  together  with  a 
transformation  for  B0.  A  suitable  transformation  may  be  defined  as 


z  rr1/2 

z  —  — ,  f  = - ,  r  =  1  -f  T]t,  B0  =  rl/7B0 

ar  a 


which,  when  substituted  into  Eq.  (35),  yields 


dB0  t  fd7B0  1  dB0  B0\ 


dr  R  \  dr2  r  dr  r2 ) 

where  R  =  ono-qa7  is  the  magnetic  Reynolds  number.  The  boundary  condition  in  Eq.  (29) 
requires  B0  to  be  continuous  at  the  interface,  which  gives 


B0  = 


at  r  =  1.  The  specific  solution  for  B0  depends  on  the  initial  condition  chosen.  In  the  present 
analysis,  two  different  cases  are  considered  for  this  condition:  an  initial  surface  current  and 
an  initial  body  current. 
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3.1  Initial  Surface  Current.  If  the  current  is  initially  located  on  the  surface  of  the 
jet,  the  appropriate  initial  condition  is  Bo  =  0  at  r  =  1.  Moreover,  since  the  boundary 
condition  given  in  Eq.  (38)  is  time  dependent,  Duhamel’s  theorem  (Carslaw  and  Jaeger 
1959)  must  be  employed  to  find  a  solution  to  Eq.  (37).  Application  of  the  theorem  requires 
the  solution  to  an  auxiliary  problem  for  Bq.  The  governing  differential  equation  for  Bq  is 
identical  to  the  equation  for  Bo,  but  the  boundary  condition  is  given  by 


K 


VqIt 

2ica 


(39) 


at  r  =  1,  where  r  is  treated  as  a  constant.  The  solution  to  this  auxiliary  problem  is  given  by 


Ba0(r,  t ,  f)  =  fr  +  £ 

J=o  yjMij) 


(40) 


where  Jo  and  J\  are  Bessel  functions,  Bq  =  “f2",  and  7_,  satisfies  the  transcendental  equation 
Ji(7j)  =  0.  This  solution  is  then  used  in  Duhamel’s  integral  to  find  a  solution  for  Bo ,  which 
gives 


where  u>j  is  defined  a s  Wj  =  7^  /(2/?)1/2 ,  and  D  is  Dawson’s  integral  (Abramowitz  and  Stegun 
1964).  In  Fie.  2,  the  dimensionless  magnetic  field  B0  is  plotted  vs.  r  at  several  times  for  a 
magnetic  Reynolds  number  of  ten,  where  B0  =  As  is  evident  from  the  figure,  the 

current  is  initially  concentrated  at  the  surface,  then  diffuses  to  the  interior  of  the  jet  as  time 
progresses.  Apparently,  for  this  magnetic  Reynolds  number  the  current  is  almost  completely 
diffused  when  r  =  2.  Also,  the  magnetic  field  at  the  surface  increases  with  time,  due  to 
the  decrease  in  the  cross  sectional  area  associated  with  the  stretching  motion  of  the  jet. 
Evidently,  this  increase  in  Bq  always  prevents  the  current  from  ever  diffusing  completely. 


3.2  Initial  Body  Current.  When  the  electrical  current  is  completely  diffused  ini¬ 
tially,  the  appropriate  initial  condition  on  Bo  is 


Voir 

2tto. 


(42) 
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Figure  2.  Bp  versus  r  at  several  times  for  R  =  10,  assuming  an  initial  surface  current. 

at  r  =  1.  Application  of  this  condition  to  Eq.  (37),  using  a  solution  procedure  similar  to  the 
initial  surface  current  case  presented  above,  gives 


Bp 


Ho  I 
2ira 


WM^j) 


(43) 


Fig.  3  shows  Bp  vs.  f  at  several  times  for  R  =  10,  where  Bp  =  •  Apparently,  for 

this  magnetic  Reynolds  number  the  electric  current  remains  almost  completely  diffused  at 
all  times.  The  positive  curvature  in  the  profiles  for  Bp  are  indicative  of  only  slight  currei  t 
diffusion  towards  the  centerline  of  the  jet.  The  slight  decrease  in  Bp  at  early  times  in  the  jet 
interior  is  an  advection  effect  associated  with  the  stretching  motion. 


4.  PERTURBATION  EQUATIONS 

To  study  the  stability  of  the  idealized  stretching  motion  presented  in  Sec.  3,  a  small 
disturbance  is  introduced  to  the  flow  and  permitted  to  develop  over  time.  The  equations 
which  govern  the  time  evolution  of  this  disturbance  are  formulated  using  linear  perturbation 
theory,  for  which  solutions  are  valid  provided  the  overall  magnitudes  of  the  disturbances 
remain  small. 
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Figure  3.  Bp  versus  r  at  several  times  for  R  =  10,  assuming  an  initial  body  current. 


The  development  of  the  perturbation  equations  is  accomplished  through  the  substitution 
of  the  variables 


Vr  =  l>rO  +  Vrl,  V*  =  V*0  +  1>*1, 


A  =  Ao  +  Aj,  p  —  po  +  Pii 
B  =  Bq  +  Bx,  BV  =  BZ  +  Bl 


n  —  ^40  +  Tj,!, 


(44) 


into  the  governing  differential  equations,  where  the  subscript  one  denotes  a  first  order  quan¬ 
tity.  Keeping  terms  which  are  of  zero-  and  first-order  only,  and  recognizing  the  zero  ordor 
solutions  also  satisfy  the  governing  equations,  the  continuity  equation  gives 


dvrl  vr,  dvzi 

~S7  +  T  +  17  =  °- 

Similarly,  the  scalar  forms  of  the  momentum  equation  yield 


(45, 
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dvT\  _  yvri  _  yr  dvr\  yz  dvr\  \  _  djh  >  dVi 
dt  2 r  2  r  r  dz  )  dr  0  dz 2 


bo 


B0 


dB_ i 


+ 


(46) 


and 


P 


yr  dvzi  y^dvzA 
2  r  dr  t  dz  ) 


^Pi  ,  /  d*vtX  ]_  dv.  1  _  d2 v2l  \ 
dz  °  \  dr 2  r  dr  dz 2  y 


dB\ 

/ip  02 


(47) 


where  substitution  of  the  zero-order  velocities  given  in  Sec.  3  has  been  performed,  and  Aj 
has  been  eliminated  using  the  relation 


r  dvz  i 
- _ — Ao. 

y  dz 


(48) 


Eq.  (48)  represents  the  first  order  contributions  to  Eq.  (20).  The  magnetic  convection- 
diffusion  equation  yields 


dB\  r}Bx  / dB0  B0\  _  yrdBx  yzdth  _  J_  f&Jh  15B,  _  B±  d2Bx\ 

dt  2 r  \  dr  r  )  rl  2r  dr  r  dz  <rpo  y  dr 2  ^  r  dr  r2  ^  dz 2  ) 

(49) 

It  is  advantageous  to  simplify  Eqs.  (45)  -  (49)  further  before  solving  them.  The  result  of 
a  previous  analysis  (Romero  1989)  has  demonstrated  a  coordinate  transformation  which 
deforms  in  accordance  with  the  stretching  motion  of  the  jet,  given  by 

z  rr1/2 

2  =  - ,  f  =  - ,  7=1+7//,  (50) 

or  a 

eliminates  the  explicit  z  dependence  in  the  perturbation  equations.  This  permits  the  ax¬ 
ial  dependence  of  the  independent  variables  to  be  Fourier  analyzed,  which  facilitates  the 
calculations.  Furthermore,  the  number  of  independent  parameters  in  the  problem  may  l  e 
minimized  if  the  perturbation  equations  are  cast  in  dimensionless  form.  Let 
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vTi  =  r/avrexp(ikz),  vzi  =  T]avzexp(ikz), 

Y - 

Aj  =  — Xexp(ikz),  pi  =  pij2a2pexp(ikz), 

V 

B\  =  ——  Bexp(ikz),  B\  =  ^—Bvexp(ikz), 
lira  2ita 

ni  =  rborbexp(ikz),  (51) 

where  k  is  the  axial  wavenumber  and  i  =  \J —  1 .  Then  the  perturbation  equations  for  the 
dimensionless  velocity,  pressure  and  magnetic  field  in  the  jet  may  be  derived  from  Eqs.  (45) 
-  (49),  and  are  given  by 


dvT  vT  ik  _ 
dr  f  +  r3/2  2 


=  0, 


dvT  vr 
T  dr  2 


At3/2 


f,dB, 
B°  df  + 


dB0  2  B0\ 
dr  +  r  ) 


(5.4) 


T 


dr 


+ 


<92t>*  1  dvz  k2 .  \ 
5r2  r  t3  2 ) 


ikAB0B, 


*(54) 


and 


dB  B_  1/2  _  £o\  .  _  t_  (d^B  }_d£_B  _  F 

dr  2r  \  dr  r  J  T  R\  dr2  f  dr  r2  r3  j 

Here  the  dimensionless  parameters  fi  and  A  are  defined  by 


(55) 


^  3 PV2a2  .  ^o/2 

K  ’  4n2pri2a4  ’ 


(56) 


and  /?  is  the  magnetic  Reynolds  number,  f)  and  A  physically  represent  the  relative  magni 
tudes  of  inertial  to  plastic  forces,  and  electromagnetic  to  inertial  forces,  respectively;  and  R 
is  a  ratio  of  magnetic  convection  to  diffusion. 


In  the  vacuum  surrounding  the  jet,  previous  analysis  (Littlefield  and  Powell  1990)  has 
demonstrated  that  the  perturbed  magnetic  field,  B",  is  identically  zero.  This  result  also 


1-1 


follows  directly  from  Ampere’s  Law,  since  the  total  current  in  the  jet  remains  constant  over 
time. 

The  boundary  conditions  applicable  at  the  interface  between  the  jet  and  the  vacuum  must 
also  be  specified,  and  may  be  derived  from  Eqs.  (25)  -  (29).  The  first  order  contributions 
to  these  conditions  are  evaluated,  keeping  in  mind  that  the  zero  order  contributions  are 
identically  satisfied  by  the  zero  order  solutions.  Furthermore,  the  first  order  contribution  of 
the  zero  order  solutions  evaluated  at  the  perturbed  boundary  must  also  be  included.  For 
example,  the  first  order  contributions  to  the  condition  F  =  0  would  be  determined  as 


E,+ 


dF0 


dr 


n>i  =  0, 


=r60 


(51) 


evaluated  at  r  =  r^o-  Consequently,  the  boundary  condition  on  B,  which  is  derived  from 
Eq.  (29),  is  given  by 


B  +  ^r,/2  +  fb  =  0,  (58) 

evaluated  at  f  =  1.  Here  the  appropriate  dimensionless  variables  have  been  used,  and 
the  axial  dependence  has  been  Fourier  analyzed.  The  remaining  boundary  conditions  are 
determined  from  Eqs.  (25)  -  (28),  and  are  given  by 


-  (— 
\4r3 


+  A 


r*/2  + 


dBo 


Or 


n  - 


r3/2 

TT 


i  dvz 


(dvr  ik  \ 

(59) 

V2  dr  +  7-3/2 “  °’ 

r  =  ° 

(60) 

and 


when  evaluated  at  f  —  1. 


dh 

dr 


r'l’v,  =  0, 


(61) 


A  suitable  set  of  initial  conditions  must  also  be  specified  before  the  perturbation  equations 
can  be  solved.  In  the  present  analysis,  the  initial  amplitude  of  the  perturbed  axial  velocity, 
was  assumed  constant.  This  gives 
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i>2  =  iC, 


(62) 


where  C  is  real.  This  specification  satisfies  the  conditions  at  the  axis  and,  apparently,  is 
sufficient  to  determine  the  initial  conditions  on  most  of  the  remaining  variables.  Eq.  (52) 
requires  the  perturbed  radial  velocity  to  be 


kCf 


and  Eq.  (61)  yields 


(63) 


kC 


(64) 


for  the  perturbed  radius.  The  initial  condition  for  the  perturbed  magnetic  field  is  different, 
depending  on  the  solution  specified  for  Bq.  For  the  initial  surface  current  case,  B  was 
assumed  zero,  and  for  the  initial  body  current  case,  B  was  assumed  to  be  a  linear  function 
of  r.  This  gives 


B  = 


kCf 

3 


(65) 


for  the  initial  body  current  case.  The  first-order  pressure  satisfies  a  Poisson-type  equation, 
which  is  formulated  from  the  time  derivative  of  the  continuity  equation.  The  solution  for  ft 
is  different  depending  on  the  solution  chosen  for  B.  For  the  initial  surface  current  case,  p  is 
given  by 


p  =  Cilo(kf)  —  C 


3ft  -  2 k2 
kn  ' 


(66) 


where  70  is  a  modified  Bessel  function  and 


Ci  = 


3ft -2 k2  (3 


Arft 


-  +  (-*) 


n 


C 


W 


(67) 


For  the  initial  body  current  case,  p  is  given  by 


A(4  -  Pf>)C  (2*>  -  3S1)C 
P  =  - TT  1 - r>7 - •"  6-2 Jo(kr) 


3k 


Qk 


(68) 


where 
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(69) 


C2 


k  2A(2  -  k2) 
8  ”  3k 


2k2  -  m 


Qk 


Io(k)‘ 


The  constant  C  is  arbitrary  and  was  chosen  as  C  =  f  so  that  =  1  initially.  Using  these 
initial  conditions,  Eqs.  (52)  -  (55),  together  with  the  boundary  conditions  in  Eqs.  (58)  -  (61), 
may  be  solved  for  the  five  unknowns  in  the  jet:  vr,  vz,  p ,  B,  and  r 


5.  SOLUTION  TO  PERTURBATION  EQUATIONS  AND  DISCUSSION 

Once  solutions  to  the  perturbation  equations  are  obtained,  they  may  be  used  to  identify 
disturbances  that  are  unstable.  A  particular  disturbance  is  considered  unstable  if  it  grows  in 
an  unbounded  manner  as  time  progresses.  Since  it  is  not  practical  to  monitor  the  solutions 
of  every  dependent  variable  when  a  disturbance  is  introduced,  the  growth  rate  of  a  single 
variable  is  often  used  to  determine  the  stability  criterion.  Fortunately,  the  selection  of  this 
variable  is  usually  not  very  critical,  because  it  is  often  obvious  when  a  particular  disturbance 
is  unstable.  In  the  present  analysis,  the  relative  perturbed  amplitude,  fy,,  is  monitored  to 
establish  the  stability  criterion.  If  |f(,|  grows  with  time,  the  disturbance  is  identified  as 
unstable. 

The  numerical  integration  ff  Eqs.  (52)  -  (55)  may  be  accomplished  in  a  straightforwar  1 
manner  using  standard  finite  difference  techniques.  A  radially  implicit  scheme  second  ordc  r 
in  f  and  first  order  in  r  was  employed.  Although  the  application  of  this  scheme  proved  to 
be  satisfactory  in  most  respects,  calculations  of  sufficient  accuracy  were  sometimes  difficult 
to  obtain.  Apparently,  there  were  two  major  causes  for  this  loss  in  accuracy.  First,  the  nu¬ 
merical  evaluation  of  the  infinite  series  in  Bo,  as  given  by  Eqs.  (41)  and  (43),  often  resulten 
in  large  truncation  errors,  particularly  at  early  times.  Consequently,  it  was  sometimes  nec¬ 
essary  to  approximate  Bo  by  direct  numerical  integration  of  Eq.  (37).  Second,  the  numerical 
evaluation  of  B  was  often  inaccurate  at  early  times  due  to  the  near-discontinuity  in  its  slope 
near  the  boundary.  As  a  result,  very  small  times  steps  were  sometimes  required  to  obtain 
accurate  solutions. 

The  result  of  a  calculation  for  f&,  using  the  initial  conditions  for  the  initial  surface  cui- 
rent  case,  is  depicted  in  Fig.  4.  The  relative  amplitude  is  shown  vs.  /,  where  t  =  rjt.  Th  ? 
parameters  fi  and  A  were  taken  to  be  one,  and  the  magnetic  Reynolds  number  was  taken 
as  ten.  Several  curves  are  shown  for  various  values  for  the  wavenumber  k.  which  is  related 
to  the  wavelength  of  the  perturbation  at  any  given  time  by  Xw  =  As  is  evident  from 

the  figure,  disturbances  of  all  wavelengths  appear  to  be  unstable.  The  most  unstable  dis- 
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turbance  at  early  times  has  an  initial  wavelength  approximately  equal  to  the  jet  diameter, 
corresponding  to  k  =  n.  At  later  times,  the  slopes  in  the  curves  for  the  shorter  wavelength 
disturbances  increase,  indicating  an  increase  in  growth  rate  for  these  perturbations.  Pre¬ 
sumably,  this  increase  would  eventually  result  in  larger  relative  amplitudes  for  the  initially 
shorter  wavelength  disturbances.  These  observations  are  in  qualitative  agreement  with  the 
instabilities  found  in  perfectly  conducting  jets  (Littlefield  and  Powell  1990). 


Figure  4.  Relative  amplitude  fb  versus  dimensionless  time  for  different  values  of  k  with 
fl  =  A  =  1  and  R  =  10,  assuming  an  initial  surface  current. 

If  the  initial  current  distribution  corresponds  to  a  uniform  body  current,  a  similar  result 
is  obtained  for  the  relative  amplitude.  In  Fig.  5  the  relative  amplitude  is  shown  vs.  r, 
assuming  an  initial  body  current,  using  the  same  parameters  previously  specified  in  Fig.  4. 
A  comparison  of  Fig.  4  with  Fig.  5  reveals  virtually  identical  behavior  in  the  growth  of  tv 
This  is  a  surprising  result  considering  the  major  differences  in  current  distribution  for  the 
two  cases,  as  indicated  in  Figs.  2  and  3.  Apparently,  a  change  in  the  current  distribution 
has  only  a  minor  effect  on  the  growth  rate  of  the  relative  amplitude. 

To  assess  the  overall  effect  of  electrical  conductivity  on  the  growth  rate  of  perturbations, 
the  relative  amplitude  fj  is  plotted  in  Fig.  6  vs.  t  for  three  different  magnetic  Reynolds 
numbers.  The  parameters  fl  and  A  were  taken  as  one,  and  the  wavenumber  k  was  taken  as 
7 r,  which  corresponds  to  the  most  unstable  wavelength  in  all  cases.  Evidently,  the  largest 
growth  rate  in  fj  is  observed  when  the  electrical  conductivity  is  very  large,  corresponding  to 


18 


fi  =  A  =  1  and  R  =  10,  assuming  an  initial  body  current. 


R  — ►  oo.  The  smallest  growth  rate  in  relative  amplitude  occurs  in  the  absence  of  electrical 
current,  or  when  R  =  0.  When  R  =  10,  the  range  in  ff,  falls  between  these  two  extremes, 
but  is  reasonably  close  to  the  values  for  R  — >  oo.  An  increase  in  the  electrical  conductivity, 
apparently,  results  in  only  small  increases  in  the  relative  amplitude. 

The  physical  mechanisms  responsible  for  the  hydrodynamic  and  MHD  instabilities  found 
in  jets  have  been  discussed  in  detail  elsewhere  (Powell  and  Littlefield  1990;  Littlefield  and 
Powell  1990)  for  jets  of  infinite  conductivity.  Evidently,  the  same  physical  mechanisms  are 
also  applicable  to  jets  of  finite  conductivity,  and  result  in  the  same  qualitative  conclusions 
about  the  stability.  In  particular,  the  application  of  electrical  current  to  a  perturbed  jet 
results  in  variations  in  the  magnetic  field  and  current  density,  which  both  increase  in  regions 
where  localized  necking  occurs.  The  current  and  magnetic  field  interact  to  create  Lorentz 
forces  that  always  decrease  the  stability  of  the  jet.  The  exact  location  and  magnitude  of  these 
Lorentz  forces  depends  on  many  factors,  including  the  assumption  for  electrical  conductivity. 
If  infinite  conductivity  is  assumed,  for  example,  the  Lorentz  forces  are  concentrated  at  the 
surface  of  the  jet.  For  jets  of  finite  conductivity,  however,  the  Lorentz  forces  are  distributed 
throughout  the  cross  section. 

Disturbances  in  the  current  clearly  result  in  decreases  of  the  stability  of  jets.  A  qualita 
tive  discussion  of  the  Lorentz  forces,  however,  does  not  sufficiently  explain  the  differences  in 
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t 

Figure  0.  Comparison  of  relative  amplitudes  at  different  times  when  R  —*  oo,  R  =  10  and 
R  =  0,  with  0  =  A  =  1  and  k  =  it. 


relative  amplitude  depicted  in  Fig.  6.  It  is  necessary  to  understand  the  effect  of  current  dis¬ 
tribution  on  the  magnitude  of  the  Lorentz  forces.  A  convenient  method  for  making  objective 
comparisons  is  to  consider  the  net  axial  Lorentz  force  per  unit  length  on  a  differential  disk 
of  width  dz  in  the  jet.  It  is  this  force  that  causes  the  disk  to  be  accelerated  out  of  the  neck. 
A  simple  computation  reveals  this  force  is  initially  given  by 

„  . 

F=— — - exp(ifcz)  (70) 

Z47ra 

if  the  current  is  concentrated  at  the  surface.  For  a  uniform  body  current,  the  axial  force  per 
unit  volume  is  initially  given  by 


/  = 


ik'CiipPr2 

127r2a5 


exp(i&£) 


(71) 


which,  when  integrated  over  the  area  of  the  disk,  gives  the  result  in  Eq.  (70).  At  later  times, 
initial  surface  currents  remain  on  the  surface  for  <t  — ♦  oo,  and  uniform  body  currents  remain 
relatively  uniform  for  finite  conductivity.  As  a  consequence,  the  axial  Lorentz  forces,  as  well 
as  the  relative  amplitude,  should  remain  approximately  equal  at  later  times  for  the  different 
current  distributions.  The  infinite  conductivity  case,  apparently,  results  in  only  slightl  v 
higher  growth  rates  for  fj,  because  the  Lorentz  forces  are  concentrated  on  the  surface  of 
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the  jet.  This  concentration  of  forces  slightly  increases  the  velocities  at  the  surface,  which 
increases  the  relative  amplitude. 

The  effect  of  current  distribution  on  the  Lorentz  forces  in  the  jet  becomes  readily  apparent 
upon  examination  of  the  velocity  disturbances.  In  Fig.  7,  the  perturbation  streamlines  are 
shown  for  a  wavelength  segment  of  the  jet  at  time  t  =  0.5,  where  ft  =  A  =  1,  R  =  10, 
and  k  =  3jt/2.  The  left  portion  of  the  figure  shows  the  streamlines  for  an  initial  surface 
current,  whereas  the  right  portion  corresponds  to  an  initial  body  current.  In  both  figures, 
a  stabilizing  secondary  motion  is  observed  in  the  central  portion  of  the  jet.  This  secondary 
motion  is  apparently  caused  by  restoring  plastic  forces  in  the  jet  (Littlefield  and  Powell  1990) 
and  is  present  even  in  the  absence  of  electromagnetic  effects.  The  initial  body  current  case, 
however,  displays  an  additional ‘cell  near  the  centerline  of  the  jet.  Evidently,  the  presence  of 
this  cell  is  caused  by  electromagnetic  body  forces  that  act  near  the  centerline  when  current 
is  distributed  uniformly  throughout  the  jet.  Even  though  the  perturbation  in  the  boundary 
appears  to  be  equal  in  both  cases,  the  disturbance  velocity  distributions  are  quite  different. 


Figure  7.  Perturbation  velocity  streamlines  for  a  wavelength  segment  of  the  jet  at  t  =  0.5 
with  k  =  3ff/2,  ft  =  A  =  1  and  R  =  10  for  (a)  an  initial  surface  current  and  (b) 
an  initial  body  current. 
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6.  CONCLUSIONS 


The  overall  effect  of  finite  electrical  conductivity  on  the  stability  of  uniformly  elongating 
plastic  jets  appears  to  be  mostly  quantitative.  The  largest  growth  rates  in  the  boundary 
perturbation  are  observed  for  jets  of  infinite  conductivity;  the  magnitude  in  fb  for  jets  of 
finite  conductivity  differs  only  slightly  from  these  results.  These  similarities,  however,  are  not 
totally  unexpected.  The  electrical  current  was  assumed  equal  in  every  case  and,  consequently 
the  average  magnetic  pressure  in  the  jet  is  about  the  same  for  all  conductivities.  The  local 
value  for  the  pressure,  however,  is  variable  and  depends  on  the  current  distribution  in  the 
jet,  as  well  as  the  conductivity.  Apparently,  the  small  differences  in  fb  result  from  small  local 
variations  in  the  magnetic  pressure.  This  pressure  becomes  concentrated  at  the  surface  for 
large  conductivities  which,  evidently,  results  in  slightly  larger  growth  rates  for  fb. 

Application  of  the  analysis  to  shaped-charge  jets  indicates  that  electrical  currents  signif¬ 
icantly  increase  the  instabilities  in  the  jets  and  may  reduce  the  breakup  time.  Calculations 
for  realistic  levels  of  electrical  conductivity  do  not  appreciably  alter  this  result.  The  increase 
in  instability  is  present  even  for  sufficiently  low  current  levels,  where  large  scale  changes  in 
the  density  and  temperature  of  the  jet  are  not  expected. 
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